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Abstract 

In the framework of the Ginzburg-Landau equation, the tempera- 
ture dependence of the upper critical field of small ring-like supercon- 
ductors is studied. At equilibrium small parts of the phase diagram 
show paramagnetism for width / radius ratios below 0.85. Their num- 
ber and extension increase with the size of the hole. In these regions, 
only the inner part of the ring shows a positive magnetic moment. 
The order parameter density profile appears to change, when crossing 
a first order transition line, which separates different angular momen- 
tum values, and we clarify the relationship between the localization of 
superconductivity nucleation and paramagnetism of those samples. 



Introduction 



The paramagnetic Meissner effect (PME) observed in small superconducting 
samples simply expresses the fact that the energy of those samples decreases 
upon an increasing magnetic field. A reentrant (H,T) phase diagram is one 
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of the signals that characterizes this effect. This surprising effect has been 
seen both in conventional [0J @ and in high Tc superconductors ||. 

The behavior of nanometer-scale superconducting grains which shows 
Pauli paramagnetism!! in addition to other striking properties of finite elec- 
tron systems, is outside the scope of this work which relies on the Ginzburg- 
Landau equation. 

Metastability is one of the major ingredients that have been put forward 
to explain PME in conventional superconductors ||. For example, in the 
Ginzburg-Landau approach of small size cylindrical samples at the upper 
critical field, the disbalance of the screening currents leads to a positive 
magnetic moment 0. But, in this case if the equilibrium transitions are 
only allowed, the PME disapears j?J, and many mechanisms causing non 
equilibrium flux configuration have been proposed ||. 

However PME is known to occur for minimal energy configurations of 
some mesoscopic objects: the Little-Parks || ring is the most simple ex- 
ample that shows an alternating sequence of diamagnetic and paramagnetic 
responses as the magnetic field is increased. On the theoretical side this gen- 
eral behavior is known to be shared by all loop-like samples and moreover 
that PME appears the be stable against thermic fluctuations H]. 

In the framework of the London approximation of the vortex phase of 
small superconducting disks, we have obtained magnetization curves which 
showed a PME, in good agreeement with the experimental observation 0, 
and made the analogy with the Little-Parks ring to interpret such a behavior 



of stable configurations [10 



The aim of this paper is to give the reader a simple example of a meso- 
scopic ring-like superconducting film of which the paramagnetic response of 
the stable configurations can be tuned with the size of the central hole, and 
to show at which part of the ring (inner or outer) the nucleation of super- 
conductivity is localized. 

In Section 2 the upper critical field H* 3 is calculated in the framework of 
the Ginzburg-Landau equation, and, if the size of the hole is large enough, the 
parts of the phase diagram where PME takes place are found. In Section 3 
we show the two regions of the ring with opposite supercurrents which might 
give rise to a global positive magnetic moment. In Section 4 we give in more 
details our results for a ring with a small hole with width / Radius = 0.75 
for which there is only one paramagnetic domain and we also show the order 
parameter density profile that changes drastically when crossing the first 
order line. Our conclusions are drawn in the last section. 
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2 Nucleation of superconductivity near the 
edges of a ring 



We shall follow quite closely the method of Bezryadin, Buzdin and Pannetier 
IT] in the case of a circular hole or of a disk[0, in order to obtain the phase 



diagram of a circular ring of external radius R and width d under a uniform 
magnetic field H normal to the plane of the film sample. 

In the case of a thin film of thickness e the effective screening length, 
X e ff = X 2 /e , becomes very large and we may neglect the magnetic field 
energy in the Ginzburg-Landau (GL) free energy functional: 



a 1*1 



- * 4 H 

2 4m 



-> 2e- . , 
-ihV A i * 



dV 



where a = a(T — T ). 

Moreover due to the small size of the sample, we can take for the vector 
potential A the potential of the uniform external applied field H : 



A = -(H x "r ) 



Using polar coordinates (p , ip), the equation for the normalized order 
parameter ip is: 



d 2 

—^ip + — r f 
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where x = p/R , ip = * / yf\a\ / b , ^ 2 = h 2 /4m \a\, 4>q = nhc/ e, 
and = 7r R 2 H is the flux of the magnetic field across the plain disk (with 
no hole). 

In order to obtain the upper critical field H* 3 [I3|, using the linearized 
version of the GL equation, we search for solutions of definite angular mo- 
mentum n of the form: 
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Then, in the present gauge, the super-current is tangent to the two cir- 
cular insulating boundaries of the ring, and we may take the boundary con- 
ditions of the GL equation as follows: 

d d 
— f\ x =i = , and -j-f\ x =i-d/R = . 

After reduction to Kummers's equation, the general solution can be ex- 
pressed as the linear combination of two confluent hypergeometric functions 



M and U 14 as follows: 



f n {x) = x n exp(-^) I. A M(Y, n + 1, x 2 ^-) + B U(Y, n + 1, x 2 ^] 

2 00 I 00 00 



where Y — ^ — \tt^t , and where the two constants A and B have to 

2 4 ^ <f> 

be determined using the boundary conditions. 

In practice, we have first obtained the ratio A/ Bin terms of , R 2 / £ 2 and 
n using the first condition, while the second condition has been numerically 
solved and allows us to obtain R 2 / £ 2 versus <fi for each n. Then, for a 
given flux , the critical line T c (0) is obtained by choosing the angular 
momentum value which minimizes R 2 I i 2 HTTJ. The limiting cases of either 
a plain disk or a circular hole can be reached, using the same procedure, with 
B = in the first case and with A = in the second one. The Little-Parks 



ring , can be recovered using a very small value of d/R , as shown on 
figure 1 . We see that paramagnetism starts just after the transition between 
angular momentum states J — 1 —>■ J and disapears for 0/0o > J, until 
the next transition occurs. On figure 2 we give the H* 3 line for a ring with 
d I R = 0.5. The cusps, associated to jumps of the angular momentum, are 
clearly seen and three small paramagnetic domains are located above these 
transitions. For sake of comparison the full disk H* 3 curve is shown on figure 
3 , on which no such domains are seen; such oscillations have been observed 
on micron-sized superconducting disks JT2]]. Our results give d / R ~ 0.85 



as the value under which paramagnetic domains are observed. 
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3 Magnetic moment and supercurrents 

In the state /„ of angular momentum n, the magnetic moment M = —OF / dH 
can be written as: 

M = J* 1 1 n R 2 e f u m f f n (x) 2 (n - ^- x 2 ) x dx. 
{^Ay Ji-d/R 0o 

In fact the local density of magnetic moment is related to the value of 
the supercurrent: 

H e 

J = — - tfW*) A l^l 2 

4m 2mc 

of which the components are: 

T fn{x) 2 , 2 , 

J v = in — x ) and J p = 

X 0o 

For a given flux and angular momentum, and when the ring is crossed 
from the inner edge (x — 1 — d/R) up to the outer one (x — 1), three various 
situations can be encountered acording to the sign of J v : 



The radial supercurrent density is negative everywhere on the ring for 
to > (i~d/R) 2 ■ gi ves a uniform standard diamagnetic response. 

The radial supercurrent density is positive everywhere on the ring for 
^ < n and gives a uniform paramagnetic response 

The current flows changes its sign at a point x = \J ■ Thus one 
can distinguish two different regions on the ring: the inner part which 
contibutes positively to the total magnetic moment and the outer part 
which gives a diamagnetic contribution, with currents in the two regions 
flowing in opposite directions. 



Let us comment Figures 1-3 on this point. The ideal Little-Parks ring 
of Figl (with d/R « 0) shows a sequence of uniform responses: J = is 
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diamagnetic, whereas the reentrant part of J = 1 is paramagnetic, etc... and 
the supercurrent changes its sign for all half-integer values of ^ . Fig2, which 
depicts the behavior of a thicker ring (with d/R = 0.5), shows that apart 
the standard uniform diamagnetic response of the J = state, the J = 1 
state gives a uniform paramagnetic response for small fields ^ < 1 . The 
magnetization of all the other states results of the competition between the 
inner and the outer currents. The case of the plain disk (Fig3) is interesting 
although no paragnetism is observed on the global magnetization: all the 
non-zero angular momentum states have a central paramagnetic response 
and a larger diamagnetic one located on the periphery of the disk. 



4 Paramagnetism in a disk with a small hole 



In this section, we have choosen to treat more extensively the particular case 
of a small hole, d / R — 0.75 , by solving the non-linear Ginzburg-Landau 
equation. This allows us to enter the superconducting region, and compute 
the densities of both the order parameter and the magnetization . In fact 
it can be shown that mixture of different angular momentum states can 
be neglected if only the lowest energy is searched for, so that we shall not 
consider such mixtures. 

Because the two boundary conditions restrict ip at two different points, 
namely x — 1 — d/R and x — 1, we have solved the non-linear GL equation 
in two steps. We have used standard numerical routines which allowed us to 
obtain f n (x) given the initial values / n (l) and df n /dx\ x=1 = , then / n (l) is 
found by solving df n /dx\ x=1 _ d/R = for / n (l). 

The H* 3 line for the ring with d/R = 0.75 is drawn on Fig. 4. Only 
one small paramagnetic domain shows up for fluxes 1.4 < 0/0o ^ 1-8 with 
J = 1. On the inset of Fig 4, which is a zoom of the paramagnetic region, we 
have drawn a fixed temperature path (A-B) with R 2 /C, 2 = 1.2 which crosses 
the first order transition line separating the J = and J = 1 phases. A fixed 
flux path (C-B) with 0/0o = 1-625 which lies entirely in the paramagnetic 
phase is also shown. The field dependence of the magnetization versus flux 
along the path (A-B) is given on Fig 5a, on which one sees the jump of the 
magnetic moment to a positive value when the first order line is crossed. This 
behavior should be compared with the smooth one of the magnetization along 
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path (C-B), in the paramagnetic phase at fixed field and versus temperature, 
on Fig 5b. 

We now give on Fig 6 the magnetization density at points A and B versus 
p/R , i.e. when the ring is crossed from the inner circle to the outer one. As 
already explained in the prevoius section, one clearly sees that the inner part 
is paramagnetic while the outer one is diamagnetic and that the balance of 
the two results in a small PME effect. Let us stress that, following this line 
of argument, and as said for disks in the previous section, one should observe 
for disks with very small holes a small domain with a positive local magnetic 
moment (around the hole) while there is no global PME effect. 

We discuss now the order parameter distribution across the ring. On Fig 
7, comparing these densities, we find that the order parameter is larger in 
the inner part than in the outer part , at point A i.e. in the diamagnetic 
region. On the contrary, in the paramagnetic region, at points B , the order 
parameter is 30% larger in the vicinity of the outer edge than around the 
hole. 

Figure 8 shows how the order parameter density profile changes when 
crossing the first order transition at 0/0o = 1-4 at fixed temperature R 2 /C, 2 = 
1.2. We have plotted the values of \ipj\ on the two, inner and outer, bound- 
aries of the ring for J = before the transition and for J = 1 after. 

It is worthwhile to notice that the regions of the ring where the order 
parameter is the largest, changes as the first order line is crossed, passing 
from the inner part before the transition to the outer one after. We have 
checked on the example of d / R — 0.5 that this result should be general. In 
this last case it is observed that, sweeping the temperature range of J = 1 
from one transition point (J = — > J = 1) to the next one (J = 1 — > J = 2), 
the order parameter profile evolves from more populated on the inner 
edge to more populated on the outer one. Let us note that, even though 
the superconducting parameter is maximal in this last region, a paramag- 
netic effect arises because of the large inner part of the disk in which the 
magnetization density is positive (see Fig. 6). 

5 Conclusions 

The paramagnetism of a mesoscopic object can obviously be due to any mech- 
anism leading to non equilibrium configurations. However, at equilibrium, 
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the PME depends strongly on the geometry of the sample and presumably 
can be observed in many cases much more involded than the toy model pre- 
sented here. Here we have shown that the giant vortex state of a pierced 
disk can exhibit PME in minimal energy configurations. The size of the 
hole controls both the number and extension of the paramagnetic domains. 
Increase of the field at fixed temperaure to enter those domains causes the 
order parameter to jump from inner to outer localization. The present ap- 
proach does not explain the PME seen on the surface superconducting states 
of 1:4 elongated ellipses |1| for which the PME in an asymmetric ring should 
be studied. Paramagnetic Meissner effect in a multiply-connected array of 
Josephson junctions has been reported [16] for which the diamagnetic current 
flows on the exterior plaquettes wheras the paramagnetic current flows in the 
inside of the sample; this is very similar to that we have described here. 
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Figure captions 



Figure 1. The phase diagram of a very thin ring with d/R = 0.01, versus 
normalized temperaure R 2 / £ 2 ~ 1 — T / T c and flux <fi / 0o i n flux quantum 
units. The surface superconducting critical line Hc 3 has been drawn for the 
first three angular momentum values. 

Figure 2. The upper critical field Hc 3 for a ring with d/R = 0.5. The 
first order lines which separate domains of different angular momentum are 
drawn as thin full vertical lines, whereas vertical dashed lines mark the ex- 
pected boundary of each paramagnetic region. The present study is valid to 
the right hand side of the Hc2 line which has been shown for comparison ( 
for the bulk material) . 

Figure 3. The upper critical field Hc^ for full disk with d/R = 1. 

Figure 4. The phase diagram of a disk with a small hole d/R = 0.75, in 
the vicinity of the upper critical fileld if C3 . The inset plot shows a zoom of 
the region 1.2 < 0/0 o < 1.8 in which a paramagnetic effect is expected. 

Figure 5a. The magnetization at fixed normalized temperature R 2 / £ 2 = 
1.2 versus flux along the path A-B ( See the inset graph of Fig. 4 ). The 
magnetization jumps to small positive values after the transition. 

Figure 5b. The magnetization in the paramagnetic region at fixed flux 
0/0o ^ 1-6 versus temperature, along the path C-B ( See the inset graph 
of Fig. 4 ). 

Figure 6. The magnetization distribution across the ring, at point A 
(full line) and at point B (dashed line) ( for A and B see the inset graph of 
Fig. 4 ). 

Figure 7. Order parameter density at point A, before the transition , 
and at point B, in the paramagnetic region, versus the radial variable 1 — 
d/R < p/R < 1 . 

Figure 8. The order parameter density profile oscillation when crossing 
the first order transition at 0/0o — 1-4 for a ring with d/R = 0.75. The 
inner edge is more superconducting than the outer one before the transition, 
whereas the reverse situation occurs in the paramagnetic phase. 
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